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ABSTRACT 

Theoretical  relationships  describing  characteristics  of  elastic- 
plastic  impact  by  rigid  spheres  on  flat  plates  were  derived.  The  theory 
is  based  on  simple  empirical  laws  relating  the  pressure  (n)  to  the  size 
of  the  indentation  or  contact.  The  derivation  neglected  the  stress  wave 
energy.  With  the  aid  of  conservation  and  fracture  mechanics  relations, 
equations  were  derived  for  the  important  impact  characteristics  including 
depth  of  damage,  remaining  strength  and  coefficient  of  restitution.  The 
theory  includes  the  equations  of  Evans  , and  Wiederhorn  and  Lawn  which  appear 
as  special  cases. 

The  theory  was  evaluated  by  comparison  with  experimental  data  obtained 
using  the  following  combinations  of  materials: 

Target  Impacting  spheres 

Transformation  toughened  zirconia  Tungsten  carbide 

Zinc  sulfide  Glass 

Glass  Glass 

The  responses  of  the  first  two  material  combinations  were  elastic-plastic 
while  the  third  was  elastic.  Agreement  with  the  theory  was  reasonable  for 
elastic-plastic  impacts  on  transformation  toughened  zirconia  and  zinc  sulfide 
but  for  elastic  impacts  on  glass,  experimental  difficulties  prevented  the 
comparisons  except  for  the  contact  radii,  so  that  further  investigation  is 
required. 


I.  introduction 


Presently  available  theories  of  solid-solid  impact  are  based  mainly  on 

(1)  (2  3) 

the  Hertz  theory  of  contact.  Evans  ’ ' used  the  Hertz  theory  to  calcu- 

late the  load  at  various  impact  velocities  and  combined  this  result  with 
fracture  mechanics  relations  to  obtain  predictions  of  the  depth  of  damage 

and  strength  degradation  for  impacts  of  spheres  on  flat  plates.  Similarly, 

(4) 

Timoshenko  and  Goodier  derived  an  expression  for  the  contact  time  of 
elastic  spheres  impacting  at  various  velocities.  They  assumed  that  the 
contact  time  was  large  compared  with  the  natural  vibration  frequency  of  the 
spheres  so  that  these  vibrations  could  be  neglected.  Yoh-Han  Pao^  extended 
the  application  of  the  Hertz  theory  to  impact  of  small  viscoelastic  spheres 
on  a steel  plate.  Evans ^ and  Wiederhorn  and  Lawn^  have  derived  expressions 
for  the  elastic-plastic  regime  and  impact  by  sharp  particles. 

In  this  paper,  a theory  of  elastic-plastic  impact  by  rigid  spheres 
on  flat  plates  is  derived.  The  theory  is  based  on  an  empirical  law  relating 
the  pressure  to  the  size  of  the  indentation  and  contact.  The  impact  response 
was  determined  for  three  combinations  of  ceramic  materials  and  the  results 
were  compared  with  theoretical  predictions.  Lastly,  the  effect  of  variations 
in  impact  conditions  was  evaluated  by  qualitative  methods. 

Several  definitions  are  necessary.  Various  stages  of  impact  are  illus- 
trated in  Figure  1.  At  first  contact  the  velocity  of  the  impacting  sphere 
is  the  Impact  velocity  (Vq) . As  the  sphere  moves  into  the  target  it  deceler- 
ates and  comes  to  rest  as  shown  in  Figure  1(B).  The  maximum  distance  that 
the  sphere  moves  into  the  target  is  the  maximum  penetration  (zc).  At  this 


point  the  radius  of  contact  between  the  sphere  and  the  target  Is  the  maximum 

_?*d lu»  (t"c)  and  the  force  Is  a maximum  called  the  load  (?)  . As 

elastic  recovery  occurs,  the  sphere  Is  accelerated  to  the  rebounding  velocity. 

After  elastic  recovery  an  Indentation  with  a projected  radius  ( r ^ ) remains 

in  the  target.  r Is  sllghtlv  smaller  than  the  r because  of  elastic  strain 
1 c 

at  the  periphery  and  elastic  recovery.  Pamage,  In  the  of  cracks,  mav 

have  spread  from  the  Impact  site.  The  strength  controlling  dimension  of  this 
damage  Is  called  the  depth  of  damage  (c). 
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II.  THEORY  OF  ELASTIC-PLASTIC  IMPACT 


Impact  conditions  and  characterist  ics  are  Hated  in  Table  I.  In  this 

paper  impact  characteristics  are  derived  for  perpendicular  impacts  hv  rigid 

spheres  on  flat  plates  at  subsonic  velocities.  The  stress  waves  generated 

bv  impact  are  neglected.  The  derivation  is  based  on  an  expression  for  the 

* 

pressure  of  indentation  based  on  the  following  empirical  law  : 


ri  r< 

n'  - ^’(— ) 


in  which  k'  la  a constant  for  a particular  combination  of  sphere  and  plate 
expressed  in  units  of  pressure  and  i, ' is  the  slope  of  the  curve  of  tog  n'  “ 

l 

fog  k'  + £'  ?og( — ) as  in  Mever's  analysis'1’  of  hardness  test  data,  k’ 
o 

and  ?. ' are  considered  to  be  the  resistance  to  indentation  and  an  exponent, 
perhaps  related  to  work  hardening,  respectively.  There  is  a similar  assumed 
relationship  between  the  pressure  of  contact  and  the  radius  of  contact 


c f. 

n - • t r> 

r 


This  expression  takes  account  of  both  the  elasticallv  deformed  rim  around 
the  indentation  and  the  indentation  itself  and  - and  are  characteristics 


of  both  the  elastic  and  plastic  deformation. 


*See  List  of  Symbols,  Table  II. 
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TABLE  I 

Impact  Conditions  and  Characteristics 


Impact  Conditions 


Account  for  or  neglect  stress  waves 

Blunt  projectile  (sphere  or  punch)  or  sharp  projectile  (cone  or  pyramid) 
Angle  of  impact 

Deformable  or  rigid  projectile 
Subsonic  or  transsonic  velocity 
Single  particle  or  multiple  particles 


Impact  Characteristics 


Phase  I 


Phase  II 


Impact  velocity 
Radius  of  contact 
Load 

Time  of  penetration 
Depth  of  damage 
Remaining  strength 


Coefficient  of  restitution 
Radius  of  indentation 
(load) 

Time  of  elastic  recovery 

(lateral  vent  crack  size) 
(remaining  strength) 


-l.ist  of  Svmbo’s- 


« 


* 


b 

c 


°0 

C10 


8 


H 

J{  ♦ > 


t 

e 

C 

P 

t , 
e 


Diagonal  dimension  of  contact  of  a regular  four  sided  pyramidal 
indenter  (Vickers  Indenter) 

Diagonal  of  a regular  four  sided  pyramidal  Indenter 
Depth- of  damage 

Velocitv  of  longitudinal  waves  of  infinite  wavelength  in  a 

L 

bar  made  of  the  same  material  as  the  target,  cQ * (E/p) 

Velocitv  of  longitudinal  waves  of  infinite  wavelength  in  a bar 

L 

made  of  the  same  material  as  the  projectile,  cj(1  “ 

Coefficient  of  restitution 


Young's  modulus  of 
Young's  modulus  of 

*)  *>  O 


<2C  -Y  r - 4;l<5 


the  target  material 
the  projectile  material 


Vicker's  hardness 

Imaginary  part  of  the  complex  number  {•} 

Critical  stress  intensity  factor  of  the  target 

Mass  of  the  projectile 

Load  developed  during  impact 

Radius  of  contact  at  time  t (r  ■ r(t)  dt) 

0 

Radius  of  a spherical  projectile 
Maximum  radius  of  contact 
Radius  of  Indentation 

Time  of  elastic  deformation  (penetration) 

Time  of  elastic/plastic  deformation  (penetration) 
Time  of  elastic  Recovery 
Time  of  plastic  deformation 


9 


V 

vc 

w 

w 

p 

w 

w 

( 

w 


pi 


el 


z 

z 

o 

z 

c 


a 


B 


Y 

? 

n 

n' 

K 


K' 


x 


TABLE  II  (Continued) 

Time  of  contact  between  target  and  projectile 
Instantaneous  velocity 
Impact  velocity 

Energy  at  time  t (W  * ft  W(t)  dt) 

o 

Energy  of  penetration 
Plastic  energy 
Elastic  energy 
Stress  wave  energy 

Constant  characteristic  of  the  geometry  of  a crack 

t 

Penetration  at  time  t (z  = / z(t)dt) 

o 

Length  of  a conical  or  a pyramidal  indenter 
Maximum  penetration 
Depth  of  indentation 

Half  angle  of  a conical  or  pyramidal  indenter 
Angle  of  incidence  of  a projectile  on  a plate 

, /mS 

v 1-2  a 

Dummy  variable  for  integration 

Pressure  of  contact 

Pressure  of  indentation 

Resistance  to  penetration 

Resistance  to  indentation 

- m-1/5  3/5  -6/5 

2t  M r g 
flj 

3 

P c 

o 
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TABLE  II  (Continued) 


/v  °1CCL 

X(1  + ~~±) 
pc 


X (1  + 


Zb  E 
Xi, 
PE 


Exponent  characteristic  of  the  variation  of  the  pressure  of  contact 
with  increasing  load 

Exponent  related  to  work  hardening 

Density  of  the  target 

Density  of  the  projectile 

Poisson's  ratio  of  the  target 

Poisson's  ratio  of  the  projectile 

Fracture  stress 


(1 . 068)  ( i+o)  (^~)  (g) 


l-o  N ,16,6/5  , f°  XirlJ^JLlL.) 

1 F (O  ' 

o 


Geometrical  "constant"  characteristic  of  the  contact 


2 

- * 7T  r 

o 


.2 

t3  Px 


(€+4) . 

K 


£+2 


♦ 
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As  in  Meyer's  law  the  pressure  of  contact  is  assumed  to  equal  the  load 

fc 

divided  by  the  projected  contact  area  leading  to 


P » k-  n r2  (-^-)C+2 
o r 

o 


(3) 


Experimental  evidence  of  the  applicability  of  equation  (1)  to  ceramics 

(3) 

has  been  presented  by  Evans  and  Wilshaw  for  static  indentations  in  zinc 
sulfide  and  by  Kirchner  et  al^*^  for  transformation  toughened  zirconia 
(TTZrO^) . The  data  at  the  larger  contact  radii  are  more  reliable  than  those 
at  smaller  contact  radii  so  the  slopes  were  based  mainly  on  the  more  reliable 
part  of  the  data. 

Assuming  that  the  stress  wave  energy  is  negligible,  the  work  done  by 
the  impact  load  during  penetration  is  equal  to  the  incoming  kinetic  energy, 
so  that 


MV  dV  = P(z)  dz 


or 


>>  < - /;  * * ‘ 
0 o 


(4) 


which  when  integrated  yields 


, K TT  r r r . 

•>  mo  - -izr-  <r-> 

o 


(5) 


Rearranging  (5)  yields  a relation  for  radius  of  contact  in  terms  of  impact 
velocity 


We  are  mainly  interested  in  the  maximum  load.  However,  as  the  load 

Increases  after  initial  contact  p is  related  to  the  instantaneous  radius 

£±2 

of  contact  (r)  by  P = < tt  r2  (~)’+2  * < tt  r2  (— ) ^ , 

J o r o r 

o o 


12 


r [H  V2j^ 

o _ J o 
k n r 

o 


which  when  substituted  in  (3)  yields  an  expression  for  the  load 


C+2  f,+2 

„ 2 fl  M(£+4)  ,f,+4  f+4 

P » tc  tt  r [S  — “ — £]  V 

o J o 


Then,  by  substituting  Q for  the  terms  not  involving  the  velocity,  equation 
(7)  can  be  shortened  to 


2(r,+2) 
, (C+4) 


P-S1V 


The  following  fracture  mechanics  relation' 


describing  the  variation  of  the  depth  of  damage  with  the  load  for  solid- 
solid  contacts  in  which  )(  Is  a "constant"  for  a particular  combination  of 
plate  and  lndenter,  accounting  for  effects  related  to  indenter  shape,  inter- 
face friction,  etc.,  was  substituted  in  (8)  to  obtain  the  following  expression 
for  the  depth  of  damage 

c - (fiV/3  V3  (10) 

Kxc 


Assuming  a uniform  tensile  stress  is  applied  during  subsequent  use  of  the 
damaged  member  so  that  the  fracture  stress  is  described  bv 
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the  expression  for  the  remaining  strength  after  Impact  damage  waa  obtained 
bv  using  equation  (11)  to  substitute  for  c In  equation  (10)  with  the  following 
result 


2(53-2) 

Y<*n)l/3  vj  5+4 


Then,  using  the  principle  that  the  impulse  equals  the  change  in  momentum, 
the  time  of  penetration  was  determined  using 


M M dV 

P “ Q v2  (5«  ' 


which  after  integration  yields 


Note  that  the  time  of  penetration  is  zero  at  infinite  velocity  so  that 
the  constant  of  Integration  is  zero. 

Now  considering  the  second  phase  of  the  impact  process  in  which 
elastic  recovery  occurs,  expressions  for  the  coefficient  of  restitution, 
the  rebounding  velocity  and  the  time  for  elastic  recovery  can  be  derived. 
Starting  *from  equation  (1)  we  obtain 

2 r-i  f'+2 

P ■*  k 1 it  r ( V’  (15) 

o r 

o 

which  is  similar  to  equation  (3) . The  energy  of  penetration  to  the  radius  r 
is,  from  equation  (5) 


r0  . r 
r„+ 4 (r  ’ 

O 
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and  the  total  energy  of  penetration  Is 


> n r r f , 
W O.  ( +4 

P C+4  V* 

o 


while  the  energy  of  plastic  deformation  is 


•cur'  r.  ... 

W . 9. 

V \+4  'r  ' 
o 


and  the  difference  between  equation  (17)  and  (18)  is  the  elastic  energy 


k it  r r 


e1’  ’ +4 


2 [<V+‘  - ( t-i)t+4l 


To  eliminate  r^,  equation  (19)  was  written  as 


*•»  - 


and  an  expression  for  r{  was  obtained  combining  (1)  and  (15) 


ri  i rr  F-F'  F'~+2 

~ - KpO  (~ ) ’ " ] ' 2 

c o 


yielding 


. 3 £+4 

« , r0  /cvf.+A,.  f , s s,TcJ.-r.', >'+2\ 

V TiT^  {1  - 1 } 

o o 


Using  the  fact  that 


W - W - W . 

e>’  p p<\ 


and  that  the  coefficient  of  restitution  squared  is  equal  to  the  ratio  of 
the  elastic  energy  to  the  impact  energy  so  that 
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, w , w , 

2 e'  . p 
e “ ■ 1 - 

V W 

P P 

and  noting  that  the  first  term  in  the  right  of  equation  (20)  is  the  impact 

•) 

energy  (equation  17),  e~  can  be  obtained  bv  dividing  W in  equation  (20) 

e 1 

bv  17  and  substituting  for  r using  equation  (b)  to  obtain 
P c 


e - 1 1 - ( 


>4  r-r/  2(7 -r) 

i HSS±4J_,r«  v 

< ...  » O 

’ . II  •- 


The  above  equation  accounts  for  extreme  variations  in  material  proper- 
ties. For  elastic  materials,  there  is  no  indentation.  In  that  case  r^ 
tends  toward  zero  so  that  * 1 tends  toward  infinity  and  e - 1 as  would  be 
expected  when  stress  waves  are  neglected.  For  perfectly  plastic  materials 
c and  h ’ are  equal  so  the  first  term  to  the  right  of  the  minus  sign  is  one. 
Also,  t,  and  •’  1 are  equal  so  that  the  exponents  of  the  remaining  terms  are 
zero  so  that  these  terms  are  equal  to  one.  Therefore,  for  the  perfectly 

plastic  case  e » 0,  as  would  be  expected  when  there  Is  no  elastic  recovery. 

If  one  assumes  that  the  plastic  deformation  occurs  before  the  elastic 

deformation,  rather  than  the  reverse,  one  can  determine  the  time  of  plastic 

deformation  bv  integrating  equa t Ion  ( 1 3)  over  the  range  of  velocities  from 
the  impact  velocitv  to  the  rebounding  velocity  with  the  following  result 


1 7 +4 , M 


S KeV  > - V 5+4  | 

o o 


(22a) 


Then  the  time  for  elastic  recovery  is  given  bv 


cp  ” *pe 


lf-1  g V ~U  - . 


(22b) 
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The  total  time  of  contact 


v 


+ 2 t „ 1h 
• ? 


t 

c 


I £+4  i M ">> 4 ~\+4 

'T*1  iTVo  13  - e 


(23) 


The  equations  derived  above  are  summarized  in  Table  III. 

The  preceding  analvais  has  yielded  expressions  for  all  of  the  listed 

impact  characteristics  in  terms  of  impact  velocity.  The  relationship  of 

this  analysis  to  the  Hertz  theorv  as  extended  hv  Timoshenko  and  P.oodier  and 

to  the  analyses  of  Evans,  and  Wiederhorn  and  Lawn  can  bo  understood  by 

comparing  the  exponent  of  the  velocity  in  equation  (23)  and  the  exponent 

of  ~ in  equation  (4)  with  the  exponents  used  in  the  other  investigations. 

The  Hertzian  case^**’^  Is  the  same  as  the  present  analyses  if  £ ■ 1 and  k 

is  considered  to  represent  the  expression  involving  the  elastic  constants  in 

the  Hertz  theory.  The  present  analysis  reduces  to  that  of  Wiederhorn  and 

Lawn^’  for  damage  bv  sharp  indenters  when  f.  ■ 2 and  to  that  of  Evans^^  when 

£ is  very  large.  Furthermore  the  a«vrraret  rv  of  the  load  history  as  predicted 

(12) 

by  Evans  in  an  elastic-plastic  impact  is  described  bv  equation  (21)  and 
equation  (14). 


I • 


Table  HI.  Equations  for  laspact  Characteristics 
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III.  EXPERIMENTAL  EVALUATION 

The  theoretical  predictions  were  compared  with  experiments  using  three 

* 

material  combinations;  t ransformat ion  toughened  zirconla  (TT-7rO,,)  plates 

impacted  bv  tungsten  carbide  spheres  (radius  0.795  mm),  zinc  sulfide  plates 

*** 

impacted  bv  glass  spheres  (radius  1.5  mml  and  glass  plates  impacted 
bv  glass  spheres  (radius  1.5  mml . TT-TrO.,  and  zinc  sulfide  show  elastic- 
plastic  responses  whereas  the  response  of  the  glass  is  elastic,  providing 
a basis  for  comparison  of  these  two  cases. 

Impact  on  TT-TrO.,  plates 


a* 


The  contact  parameters  for  the  tungsten  carbide  spheres  against  TT-2r0, 
plates  were  determined  bv  "static"  loading  with  the  following  results 

k.  - 23.5  ONm”'  i - 1.06 
v’  - 22  ONm'2  - 0.52 

so  that  substitution  in  equation  (211  yields 

2 -2  0.43 

e - 9.37  . 10  V 

The  Impact  damage  was  characterized  previously.^  ^ The  above  equation  is 
compared  with  the  experimental  data  bv  Kirchner  et  nl^1  in  Figure  2-  There 
is  a reasonable  agreement  between  the  two  curves  but  the  theorv  predicts 
higher  values  of  e than  were  observed  bv  experiment.  As  indicated  in  the 


TT-7.r0,  plates  in  the  as-fired  condition  were  supplied  bv  R.  C.  Oat  vie 
of  CSIRO. 

** 


*** 


7.nS  supplied  bv  Raytheon  Co.,  Waltham,  Mass. 
Plate  glass  from  Centre  Class,  State  College,  Pa. 


J 


20 


E 


derivation,  the  theory  does  not  account  for  stress  wave  energy  losses. 
Therefore,  one  would  expect  the  theoretical  rebounding  velocities  to  be 
higher  than  those  measured  experimentally  which  is  consistent  with  the 
observations. 

Comparisons  were  also  made  for  the  depth  of  damage  and  remaining 
strength  as  shown  in  Figures  3 and  4.  These  comparisons  depend  on  the 
values  of  x»  T and  all  of  which  are  subject  to  some  uncertainty.  Des- 
pite this  fact,  there  is  reasonable  agreement  between  theory  and  experiment. 

It  should  be  noted  that,  as  indicated  in  Figure  9,  the  theory  does  not 
provide  for  the  threshold  velocity  necessary  to  induce  the  radial  cracks 
but  does  provide  a good  estimate  of  the  slope  of  the  crack  length  vs.  impact 
velocity  curve.  The  sudden  change  of  slope  of  the  experimental  curve 
above  70m’s*  in  Figure  3 is  probably  a real  effect  because  it  was  also 
observed  in  the  under-aged  zirconia  at  somewhat  higher  velocities.  At  present, 
no  explanation  of  this  anomaly  is  available. 

There  is  good  agreement  between  the  theoretical  and  experimental  values 
of  remaining  strength  of  as  fired  TT-ZrO^  (Figure  4)  . Based  on  higher 
values,  the  theoretical  curves  for  under-aged  and  peak-aged  TT-ZrO^  would 
be  expected  to  yield  higher  values  of  remaining  strength  than  the  as-fired 
material.  However,  < and  £ have  not  yet  been  determined  for  these  materials 
so  the  theoretical  curves  are  not  available. 

As  mentioned  above,  neglect  of  the  stress  wave  energy  in  the  theory 


Li 


leads  to  higher  values  of  e,  but  it  also  leads  to  predictions  of  the  load  . 
depth  of  damage  and  strength  degradation  that  are  too  high.  The  comparisons 
in  Figures  3 and  4 are  consistent  with  these  expectations. 


Threshold  for  radial  crack  formation 


20  40  60  80  100 


Impact  Velocity -m^s”1 

Fig.  3 Maximum  radial  crack  length  vs  impact 
velocity  for  TT-Zr02  (as  fired)  impacted 
by  tungsten  carbide  sphere. 
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Rem 


Experimental  curve  TT-Zr02 
X\ (under-aged)  (KIC  * 9.7  MPam,/2) 
\ 


20  40  60  80  100 


Impact  velocity  -ms 

Fig. 4.  Remaining  strength  vs.  impact  velocity  for  TT-ZrO 
as  fired  and  under-aged,  impacted  by  tungsten 
carbide  sphere. 


Impact  on  zinc  sulfide  plates 


The  contact  parameters  for  the  glass  spheres  against  zinc  sulfide  plates 
were  determined  bv  static  loading  with  the  following  results: 


< = 5. 

5 

ON 

-2 

m 

£ = 0.78 

K'  = 5. 

5 

ON 

-2 

m 

5'  = 0.27 

A series  of  indentations  at  loads  of  222,  445,  668  and  891  N are  shown  in 
Figure  5.  The  figure  shows  the  increase  in  indentation  radius  and  radial 
crack  length  with  increasing  load.  The  figure  also  shows  the  absence  of 
lateral  vent  cracks  at  low  loads  (no  shaded  area  between  radial  cracks)  and 
the  increasing  lateral  vent  crack  length  at  higher  loads. 

Zinc  sulfide  plates  were  impacted  by  glass  spheres  at  various  velocities. 
Indentations,  Hertzian  cone  cracks,  radial  cracks  and  lateral  vent  cracks 
were  observed.  The  damage  in  a specimen  impacted  at  68.7  m-s ^ is  shown  in 
Figure  6.  The  figure  shows  the  indentation  surrounded  by  a Hertzian  cone 
crack  which  was  partially  exposed  by  chipping  of  the  material  between  the 
cone  crack,  a lateral  vent  crack  and  the  surface.  The  radial  cracks  which 
are  not  clearly  shown  in  the  photograph  extend  considerably  beyond  the  chipped 
region. 

By  substituting  K , k',  £ and  in  equation  (21) 


[1  - 4.84  • 10 


-2 


V0-45]15 


This  equation  is  compared  with  experimental  impact  data  in  Figure  7.  Again, 
as  in  the  case  of  TT-ZrO^,  the  experimental  curve  falls  below  the  theoretical 
curve  at  least  partly  because  the  stress  wave  energy  losses  are  neglected 
in  the  theory. 

The  radial  cracks  were  assumed  to  be  the  strength  controlling  flaws. 

The  maximum  radial  crack  length  was  taken,  in  each  case,  as  the  distance  from 


P = 22  2 N 


P = 445 N 


Fig.  5.  Indentations  in  zinc  sulfide  caused  by  a spherical 

glass  indenter  ( 1.5mm  radius)  at  various  loads  ( I4X) 


P = 668  N 


P = 891  N 


Fig.  7 Theoretical  and  experimental  coefficients  of  restitution  vs 
impact  velocity  for  ZnS  plate  impacted  by  glass  spheres 
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the  edge  of  the  Indentation  to  the  tip  of  the  crack.  The  results  for  the 
static  tests,  plotted  in  accordance  with  equation  (8)  are  given  in  Figure  8. 
The  slope  of  the  curve  is  K /x  ao  that,  taking  K - 1 MPa  m\  X was 

Iv  1L 

_3 

found  to  be  54.7  • 10  . Using  this  value  of  the  theoretical  depth  of 

damage  was  calculated  and  compared  with  the  experimental  values  from  the 

impact  experiments  as  shown  in  Figure  9.  Again,  the  theory  yields  reasonable 

predictions.  The  depth  of  damage  predicted  by  the  theory  is  too  small 

which  is  the  opposite  of  the  result  expected  based  on  the  assumption  that 

the  stress  wave  energy  can  be  neglected,  so  that  another  explanation  is 

3/2 

needed.  Figure  8 shows  that  c is  increasing  faster  with  load  than 
expected  based  on  equation  (9).  Therefore,  further  investigation  of  the 
fracture  mechanics  relations  appears  to  be  necessary.  In  this  case  the 
remaining  strengths  of  the  zinc  sulfide  specimens  were  not  determined  because 
an  insufficient  number  of  specimens  was  available. 


Impact  Velocity-m*s_l 

Fig.  9 Depth  of  damage  vs  impact  velocity  for  impact 
of  glass  sphere  (3  mm  diameter)  on  ZnS  plate. 
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Impact  on  g lass  pint e s 

The  contact  parameters  for  glass  spheres  against  glass  plates  were 
determined  bv  static  loading  with  the  following  results: 

_ 2 

k - 46.9  CINm  £ - 2.04 

Class  spheres  do  not  cause  permanent  lndentat Ions  In  the  glass  plates  so 
' ' and  were  not  determined.  It  should  be  pointed  out  that,  whereas 
the  theory  assumes  a rigid  sphere,  the  contact  between  glass  spheres  and 
plates  is  an  elastic  contact  between  materials  having  approximately  equal 
values  of  Young's  modulus.  Therefore,  the  elastic  deformations  in  the  tvo 
bodies  are  approximately  equal  leading  to  a greater  contact  tadius,  less 
penetration  and  lower  pressures  than  one  would  expect  for  a rigid  sphere. 
Despite  their  apparent  simplicity,  these  contacts  between  glass  spheres  and 
glass  plates  are  far  from  ideal  in  terms  of  the  theory. 

In  the  previous  cases  the  evaluation  of  the  theory  was  introduced  bv 
comparing  theoretical  and  experimental  values  of  the  coefficient  of  resti- 
tution e.  In  the  case  of  glass  there  is  no  indentation  and  stress  waves 
are  neglected  so  the  theory  predicts  e ■*  1.  However,  e for  impacts  of  glass 
spheres  on  glass  plates  is  - 0.94  at  20  m-s  * and  decreases  verv  slowlv 
with  increasing  velocity  until  crushing  occurs  at  velocities  N 100  m-s  * . *■ ' 

The  theoretical  contact  radii  are  compared  in  Figure  10  with  another 

theoretical  curve  based  on  the  Hertz  theory  and  with  experimental  data.  The 

experimental  data  fall  near  the  theoretical  curves  but  appear  to  have  a lower 

slope.  The  effect  of  neglecting  the  stress  wave  energy  would  be  expected 

to  have  a weak  effect  on  r because  r varies  as  the  cube  root  of  the  load. 

c c 

However,  the  stress  wave  energy  increases  strongly  with  velocity  so  that  this 
effect  may  account  for  at  least  part  of  the  lower  slope  of  the  experimental 
data . 
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Impact  Velocity -m-sH 


Fiq.  10  Radius  of  contact  vs  impoct  velocity,  theoretical 

curve  and  data  points  from  impact  glass  onglossi,4l 
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The  variation  of  the  load  with  the  impact  velocity  was  compared  with 
theoretical  predictions  based  on  the  present  theory  and  the  Hertz  theorv. 

The  agreement  la  reasonable  but  the  experimental  data  are  very  scattered 

3 

because  the  loads  were  calculated  from  the  Hertz  theorv  in  which  the  P «■  r 

c 

so  that  errors  in  measurement  of  r have  a large  effect  on  P.  For  this  pur- 
pose a direct  method  of  measuring  the  load  is  highly  desirable. 

The  variation  of  the  depth  of  damage  with  impact  velocity  was  investi- 
gated. At  low  Impact  velocities  the  contact  radii  and  the  depths  of  damage 
were  of  approximately  equal  magnitude.  Equation  (9)  for  the  depth  of 
damage  is  based  on  the  assumption  that  the  crack  is  large  compared  with  the 
contact  radius  so  that  the  details  of  the  load  application  do  not  affect  the 
stress  distribution  near  the  crack  tip.  The  data  show  that  the  cracks  do  not 

become  sufficiently  large  until  impact  velocities  in  excess  of  60  m • s * are 
(14) 

reached.  At  higher  velocities,  a plot  of  Cog  c vs.  Cog  appears  to 

be  approximately  a straight  line  but  the  slope  is  much  greater  than  that 

expected  based  on  the  present  theorv  or  the  Hertz  theorv.  The  range  of 

available  data  at  these  higher  velocities  is  verv  limited  because  crushing 

occurs  with  a decrease  in  c at  velocities  over  about  115  m-s  V 

The  above  results  show  that  except  for  the  r vs.  V relationship  sub- 

c o r 

stantial  experimental  difficulties  have  so  far  prevented  a satisfactory 
comparison  of  the  present  theorv  with  experiment  for  the  case  of  impacts  of 
glass  spheres  on  glass  plates. 


IV.  DISCUSSION 


The  following  discussion  deals  mainly  with  qualitative  evaluation  of 
the  effect  of  variations  In  the  contact  parameters  and  impact  conditions. 

Slgnif  lcance  of  1 1 £ and  . 

The  significance  of  K,  k ' , £ and  £'  varies  somewhat  depending  on  the 

nature  of  the  deformation  at  the  impact  site.  The  definition  of  k (equation 

2)  indicates  that  v Is  the  pressure  existing  when  the  contact  radius  equals 

the  sphere  radius.  If  the  deformation  is  elastic,  k accounts  for  the  role 

of  the  elastic  constants  in  the  deformation.  For  example,  if  a sphere  is 

pressed  into  a plate  having  the  same  elastic  constants  the  contact  radius 

( 4)  ) 

based  on  the  Hertz  theorv  is 


r 

c 


P r , 


(24) 


Substituting  in  equation  (2)  assuming  P 


2 

it  r 


c 


11  yields 


U.3tt; 


r 


Thus,  if  £ - 1 as  one  expects  for  the  Hertzian  case,  then 


(25) 


* " 1737  <26) 

The  experimental  results,  assembled  in  Table  IV,  indicate  somewhat  higher 
values  of  *,  and  £ for  the  elastic  contact  (glass-glass)  than  one  would  expect 
based  on  the  Hertzian  analysis. 
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For  the  cases  in  which  there  is  plastic  deformation  of  the  target  or 

projectile,  < and  <’  are  related  in  some  wav  to  the  hardness  of  the  softer 

of  the  contacting  materials.  F and  £'  take  care  of  the  increase  in  the  rate 

of  increase  of  the  pressure  as  the  r/r  ratio  increases.  In  materials  that 

o 

deform  plastically,  this  variation  may  be  caused,  in  part,  by  work  hardening. 

The  constants  K,  £ and  F'  are  difficult  to  measure.  The  range  of 
pressures  covered  by  the  measurements  is  much  smaller  than  desired.  At  low 
pressures  the  scatter  in  the  results  is  large  because  of  experimental  errors. 
At  higher  pressures  the  measurements  are  likely  to  be  terminated  by  crushing 
of  the  sphere.  Greater  weight  should  be  given  to  the  measurements  at  higher 
pressures  in  determining  the  slopes  (F  and  F').  However,  at  best,  the  present 
techniques  yield  only  rough  estimates. 

Stress  Wave  Energy 

Stress  wave  energy  can  be  expected  to  have  a significant  role  in  impact 
phenomena.  In  elastic  impacts  (for  example,  glass-glass)  the  stress  wave 
energy  is  the  dominant  loss  mechanism  and  for  firmly  supported  targets  it 
is  the  principal  factor  determining  the  coefficient  of  restitution^^ . 
Hunter^^  gave  the  following  expression  for  the  stress  wave  energy  in  a 
target  as  a function  of  velocity. 

2 * ^ / S 

W ■=  *<  MV  X V ' (27) 

s 

, T -1/5  3/5  -6/5 

2 t M r g 

where  X - y ■ -°- — — — (28) 

0Co 


36 


i i 


in  which  0 is  the  Poisson's  ratio  and  the  lack  of  a subscript  signifies  the 
target  and  the  subscript  1 signifies  the  projectile  and 


T 


(1.068)5(1+o)( 


l-o  . h ,16,  6/5 
l-2o'  45; 


(^2-l)^?dCi 

FJO  1 
o 


(30) 


where  J indicates  the  imaginary  part  of  the  complex  variable.  The  details 
of  the  derivations  are  given  by  Hunter^^  . However  in  our  case  an  additional 
term  should  be  added  to  the  equation  to  account  for  the  stress  wave  energy 
in  the  projectile. 

For  a rigid  sphere,  a derivation  similar  to  that  for  a half  space  can 

be  done.  According  to  Love^*^,  a half  space  and  a sphere  are  one  and  the 

same  thing  mathematically  speaking.  They  are  geometrically  inverse  of 

each  other.  Indeed,  taking  the  derivation  of  Miller  and  Pursev^^’^^  for 

a half  space  and  modifying  it  by  applying  the  boundary  conditions  for  a 

sphere,  introduces  new  terms  in  the  differential  equations  (37-40)  in 

reference  14.  These  new  terms,  due  to  application  of  Hankel’s  transform  to 

(20-23) 

the  equations  of  elasticity  , are  dependent  upon  the  displacement  and 

distortion  of  the  surface  of  the  sphere,  as  well  as  their  derivatives  in 
regard  to  the  radial  direction  of  cylindrical  coordinates  centered  at  the 
impact  site  with  an  axis  perpendicular  to  the  surface.  With  the  assumption 
of  a rigid  sphere,  these  additional  terms  will  vanish  and  the  solutions 
given  by  Miller  and  Pursey^^’^^  and  by  Hunter  for  a half  space  will  be 
valid  for  a projectile.  Hence,  the  total  stress  wave  energy  absorbed  during 
impact  can  be  written  as 

W * X (%  MV*5)  V3/5 
s o 


(31) 


37 


in  which 


^ o 

X = X (1  + —2-) 


pc 

o 

with  the  impedance  ratio,  ~ ~ , taking  account  of  the  distribution  of 

i o. 


energy  between  the  target  and  the  projectile.  This  stress  wave  energy 
will  enter  on  the  right  side  of  equation  (5). 


h MV2  = - - • -Q  (—  + h MV2  X V3/5 

o £+4  r oo 

o 


In  the  relationships  previously  derived  (Table  II)  the  term  V will  thus 

o 

3/5 

be  replaced  by  V (1-XV  ) . Because  of  absorption  of  energy  by  the 


stress  waves,  the  radius  of  contact,  load  and  depth  of  damage  will  be  smaller 


than  indicated  by  equations  (6,8,10)  while  the  remaining  strength  and  contact 


time  are  greater  than  expected  based  on  equations  (12)  and  (23).  The 


coefficient  of  restitution  will  be  smaller  than  expected. 


Sharp  Indenter 


For  a sharp  indenter,  a similar  derivation  can  be  done  in  which 


n - *(f-)C 

O 


n = K(f)f’ 


where  zQ  is  the  height  of  the  indenter  for  a conical  indenter  and  a and  b 


are  the  diagonals  of  a pyramid  at  the  contact  and  at  the  base  as  shown  in 


Figure  11.  Thus,  the  equations  in  Table  III  can  be  adapted  for  the 


appropriate  sharp  indenter  by  substituting  the  symbols  as  shown  below: 


* (9) 

This  case  is  similar  to  the  Ludvick  or  Vicker's  hardness  . 
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Sphere  Cone 

Contact  Parameter  r r 

Normalizing  Parameter  r z 

o o 


Pyramid 

a 

b 


Sharp  lndenters  are  expected  to  increase  the  depth  of  damage 

that  of  blunt  lndenters.  Hence,  t he  strength  will  be  lower. 

( 24 ) 

have  been  observed  lor  static  indentations  “ , and  grinding 

one  would  expect  similar  results  for  impact. 


compared  with 


Such  results 

(25,26)  , 

and 


Defo r ma  b 1 e particle 

If  the  impacting  particle  is  deformable,  the  peak  load  will  be  lower 
than  for  a rigid  sphere  because  part  of  the  impact  energy  is  absorbed  bv 
the  work  of  bulk  deformation  of  the  particle.  The  elastic  deformation 
will  induce  vibrations  (like  bending  waves  in  a plate)  in  the  particle  and, 
to  the  extent  that  the  particle  carries  awav  this  energy,  it  is  an  impact 
energy  loss.  The  plastic  deformation,  if  it  occurs,  represents  another 
energy  loss  mechanism.  Consequently,  assumption  of  a rigid  projectile  over- 
estimates the  maximum  load,  the  depth  of  damage  and  the  strength  degradation 


Impact  at  an  angle 

For  impacts  at  an  angle  on  brittle  materials,  one  can  assume  that  the 

impact  damage  Is  controlled  by  the  vertical  component  of  velocity.  Therefore 

compared  with  perpendicular  impact.  Impact  at  an  angle  leads  to  lower  loads, 

depth  of  damage  and  strength  degradation.  In  their  analysis  of  erosion, 

(27) 

B.  .1.  Hockev  et  al  made  similar  assumptions. 
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V.  SUMMARY  AND  CONCLUSIONS 


A comprehensive  theory  of  elastic-plastic  impact  was  derived.  The 
theory  was  evaluated  by  experiments  involving  impacts  of  tungsten  carbide 
spheres  on  TT-ZrO^,  glass  spheres  on  zinc  sulfide  and  glass  spheres  on 
glass.  Agreement  with  the  theory  was  reasonable  for  elastic-plastic  impacts 
on  TT-ZrO^  and  zinc  sulfide  but  for  elastic  impacts  on  glass  experimental 
and  theoretical  problems  were  encountered  that  require  further  investigation. 

Approaches  by  which  the  theory  can  be  extended  to  take  stress  waves, 
sharp  particles,  deformable  particles,  and  impact  at  an  angle  into  account. 


were  described, 
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Abst  ract 


Zirconia  specimens,  partially  stabilized  with  calcia  and  aged  for 
various  times  to  vary  K ,,  were  impacted  by  glass,  steel  and  tungsten 
carbide  spheres.  The  impact  damage  was  character ized . The  zirconia 
was  resistant  to  damage  by  the  softer  spheres,  glass  and  steel,  bvit  the 
tungsten  carbide  caused  visible  damage  (indentations)  at  impact  velocities 
as  low  as  9 ms  ^ . At  higher  velocities,  radial,  lateral  vent,  and  circum- 
ferential cracks  were  observed.  Indentation  and  stress  wave  energy  loss 
mechanisms  account  for  most  of  the  observed  energy  losses.  The  remaining 
strengths  after  impact  increased  somewhat  less  than  expected  with  increas- 
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1.  Introduction 

The  fracture  toughness  of  zirconia  can  be  increased  by  proper 

control  of  composition  and  heat  treatment*1*  Calcia  partially 

stabilized  zirconia  was  used  in  this  investigation.  In  this  material 

aging  at  intermediate  temperatures,  usually  in  the  range  1200-1500°C 

(4 ) 

in  the  cubic-tetragonal  phase  field  , leads  to  coarsening  of  the 
tetragonal  precipitate  in  the  much  larger  cubic  crystals.  At  room 
temperature,  this  coarsening  increases  the  susceptibility  of  the  tetra- 
gonal phase  to  transformation  to  the  monoclinic.  Tensile  stress  further 
encourages  this  transformation  because  the  monoclinic  phase  is  less  dense 
than  the  tetragonal  phase.  Therefore,  when  cracks  propagate  through  the 
material,  the  tensile  stress  concentrated  at  the  crack  tip  causes  the 
tetragonal  crystals  within  a few  micrometers  of  the  crack  to  transform 
to  monoclinic  *'~>* . This  transformation  is  believed  to  be  responsible 
for  the  high  values  of  critical  stress  intensity  factor  (K.  ) observed 

l v> 

for  these  materials.  Therefore,  the  material  is  sometimes  termed  trans- 
formation toughened  zirconia  (TT  ZrO„). 

TT  ZrO , is  unique  in  that  machining  (grinding)  causes  an  increase 
in  strength^’*1*.  This  strength  increase  is  believed  to  be  caused  by  com- 
pressive surface  stresses  formed  by  the  expansion  during  transformation 
to  the  monoclinic  phase  which  is  induced  by  the  grinding.  The  decision 
to  investigate  localized  impact  damage  in  TT  ZrO.,  was  based  on  the  high 
K values  and  the  unique  surface  effect  in  these  materials  which  were 

« Ivj 

expected  to  reduce  the  impact  damage. 


II.  Procedures 


Calcia  partially  stabilized  ZrO,  with  four  1300°C  aging  treatments 
(as-fired,  0 hours;  underaged,  30  hours;  peak-aged,  48  hours;  and  over- 
aged, 57  hours)  was  prepared  in  the  form  of  rectangular  bars,  3 x 10  x 
40  mm,  with  ground  (50  grit)  surfaces.  One  large  surface  of  each  bar 
was  polished  with  6 um  and  1 urn  diamond  compound  and  either  1/4  pm  diamond 
compound  or  1 pm  gamma  alumina  in  successive  steps  on  a paper  covered 
rotating  lap.  At  intervals  during  polishing  the  relative  amounts  of  mono- 
clinic phase  in  the  surfaces  were  monitored  by  x-ray  diffraction.  Polish- 
ing was  stopped  when  the  reduction  of  the  monoclinic  phase  no  longer 
continued.  The  other  surfaces  of  the  bars  were  left  in  the  machined 
condition.  One  reason  for  doing  this  was  the  expectation  that  the  com- 
pressive surface  stresses  due  to  machining  would  reduce  the  probability 
of  fractures  originating  at  the  edges. 

A preliminary  investigation  was  used  to  select  a suitable  sphere 
material  and  size  for  the  impact  tests.  Class,  steel  and  tungsten  carbide 
spheres  were  tested  and  tungsten  carbide  was  chosen  for  the  main  impact 
experiments.  In  these  experiments,  ZrO.,  bars  with  each  of  the  four  treat- 
ments were  impacted  on  the  polished  surfaces  at  various  velocities.  The 
impacting  and  rebounding  velocities  were  measured  by  measuring  the 
distances  traveled  by  the  spheres  between  flashes  of  a stroboscope*  as 
indicated  by  phtographs  taken  by  a camera  with  its  lens  open  during  the 
flight  of  the  sphere.  The  coefficient  of  restitution  was  calculated  as 
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the  ratio  ot  the  rebound  velocity  tv'  tin'  impact  velocity.  The  specimens 
were  coatevl  with  a thin  lavor  ot  .soot  . The  contact  raJ  li  were  measured 
t'rom  tl>e  impressions  in  the  soot  . The  Impact  damage  was  characterised. 
IT  i'.rO , is  indented  readily  bv  the  tungsten  carbide  spheres.  The  radii 
ot  the  indentations  were  measured.  In  each  specimen  the  radius  ot  the 
longest  radial  crack,  it  present,  was  measured  as  the  distance  t'rom  the 
center  ot  the  indentat  Ion  tv'  the  tip  ot  the  crack. 

I'he  flexural  strengths  ot  t tie  impacted  specimens  were  measured  by 
three  point  loading  on  a one  inch  span  to  evaluate  t he  strength  degrada- 
tion caused  by  the  impact  damage.  I'he  tract  ure  surfaces  were  examined 
bv  optical  and  scanning  electron  microscopy  and  the  depth  ot  damage 
(.depth  ot  the  disturbed  areal  under  the  impact  site  was  measured.  I'he 
critical  stress  intensity  factors  IK  ^,1  ot  the  materials  were  measured 
by  the  single  edge  notched  beam  f SKNBl  method v<  and  estimated  from 
the  fracture  stresses  and  t lie  depths  ot  damage  or  radii  of  radial  cracks 
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111.  Results  and  Discussion 


Strengths  of  the  mat er Iais  before  impact 


A preliminary  aging  curve  was  available  for  specimens  of  this 
particular  composition  aged  at  1300°C  (Figure  1).  These  preliminary 
strengths  were  determined  using  bars  (1/8  x 1/8  x 1-1/2  in)  with  150  grit 
diamond  ground  surfaces  from  which  40  pm  was  removed  on  the  final  pass. 

The  strengths  were  measured  in  tour-point  loading  (1.125  in  outer  span, 

0.5  in  inner  span)  at  a loading  rate  that  resulted  in  failure  in  1-2  min. 
This  curve  was  used  to  select  the  aging  times. 

One  of  the  large  surfaces  of  each  specimen  was  polished  to  provide 
a smooth  surface  for  observation  of  the  impact  damage  and  to  remove  the 
compressive  surface  stresses  caused  by  grinding  so  that  these  stresses 
would  not  reduce  the  extent  of  damage  as  has  been  shown  previously  for 
other  materials  ^ . The  final  step  in  polishing  the  first  three  spec- 
imens of  each  aging  treatment  used  1 pm  gamma  alumina  in  liquid  suspension. 
This  process  reduced  the  monoclinic  phase  in  the  peak-aged  material  to 
about  32%  of  the  amount  present  in  the  original  ground  surfaces.  Sub- 
sequently, the  remainder  of  the  specimens  were  polished  by  a process  in 
which  1/4  pm  diamond  compound  was  used  as  the  last  step.  This  process 
reduced  the  monoclinic  phase  in  the  surface  to  about  10%  of  the  amount 
originally  present.  The  flexural  strengths  of  three  specimens  from  each 
of  these  aging,  grinding  and  polishing  treatments  were  measured  and  these 
results  are  including  in  Figure  1.  The  strengths  were  measured  in  three 
point  loading  on  a one  inch  span  with  the  polished  surfaces  in  tension. 


• Average  Flexural  Strength, As-Ground  Condition 
o Average  Flexural  Strength  .Polished  (Alumina) 
x Average  Flexural  Strength , Polished  (1/4 pm  Diamond) 


-As  Fired 
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Removal  ol  tin*  ground  surface  decreased  the  strength,  as  expected.  The 
decrease  was  greater  lor  specimens  polished  with  1/A  pm  diamond  than  it 
was  lot  those  polished  witli  ! tint  alumina  which  Is  consistent  with  the 
traction  ol  the  monoclinic  phase  remaining  in  the  surface  after  polishing. 
Removal  ot  the  stressed  surface  reduces  the  strength  hy  up  to  about  100 
MPa.  This  decrease  can  he  compared  with  decreases  of  about  (>()  and  120  MPa 
observed  bv  Pascoe  and  Carvle^’^  in  two  sets  of  experiments.  The  strengths 
of  the  over-aged  specimens  were  sometimes  very  low  because  of  the  presence 
ot  cracks  caused  by  excessive  formation  ol  the  monoclinic  phase. 

The  values,  determined  bv  the  SKN11  method,  were  as  follows: 

As- fired  5. HA  + 0.77  MPam1 ^ 

Under-aged  H.71  + 0.b7 

Peak-aged  l 2 . 99 

This  method  tends  to  yield  values  that  are  slightly  higher  than  those 
obtained  by  other  methods,  even  witli  more  conventional  ceramics.  In  the 
present  case,  the  compressive  stresses  induced  by  the  phase  transformation 
during  sawing  of  the  notch  can  be  expected  to  cause  additional  increases 
in  the  measured  values.  These  Increases  should  be  expected  to  Increase 
with  increasing  susceptibility  to  the  tetragonal  to  monocllnlc  phase 
transformation.  The  measured  values  of  K can  be  reduced  bv  correcting 
for  the  effect  of  the  residual  surface  stress.  If  It  Is  assumed  that  the 
residual  stress  In  the  notch  Is  approximately  equal  to  the  residual  stress 
originally  Induced  in  the  surface  bv  grinding  and  subsequently  removed  by 
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polishing,  the  ratio  of  the  strength  in  the  polished  condition  to  that 
in  the  ground  condition  (Figure  1)  can  be  used  to  correct  the  above  K. 

1 L 

values  with  the  following  results: 

As- fired  5.5  MPam^*" 

Under-aged  7.6 


Peak-aged 


10.5 


f 


These  values  will  be  compared  with  values  based  on  flaw  size  estimates 
in  a later  section. 
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Preliminary  experiments  to  select  impacting  spheres 


TT  ZrO.,  specimens  were  impacted  by  glass,  steel  and  tungsten  carbide 
spheres  to  provide  information  on  which  to  base  the  selection  of  spheres 
to  be  used  for  the  principal  experiments.  Class  spheres  (3  mm  diam)  did 
not  induce  visible  damage  in  the  surfaces  at  velocities  up  to  124  ms  ^ 
at  which  the  glass  crushed  completely.  Steel  spheres  (1.59  mm  diam) 
did  not  induce  visible  damage  at  velocities  up  to  198  ms  ^ but  the 
spheres  flattened  so  that  the  use  of  higher  velocities  seemed  unpromising 
Tungsten  carbide  spheres  (1.59  nun  diam)  Indented  the  surfaces  readily 
even  at  low  velocities  and  also  caused  crack  formation.  The  coefficients 
of  restitution  of  the  three  materials  were  consistent  with  these  observe- 
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It  was  clear  from  these  experiments  that  TT  Zrl^  was  very  resistant 
to  impacts  by  the  softer  projectiies,  glass  and  steel.  In  many  ceramics 
Hertzian  cone  cracks,  visible  as  ring  cracks  in  the  surfaces,  or  other 
visibLe  damage,  would  have  been  induced  by  these  spheres.  In  order  to 
be  able  to  observe  and  characterize  variations  in  localized  impact  damage 
in  the  available  range  of  velocities,  tungsten  carbide  spheres  were 
selected  for  the  principal  impact  experiments. 


Cl lar a e ter  i z a t ion  o f impact  d amage 


The  types  of  damage  observed  at  the  impact  sites  included  indenta- 
tions and  radial,  lateral  vent,  and  circumferential  cracks.  The  impact 
sites  for  two  impact  velocities  are  shown  in  Figure  2.  At  13.9  ms  ^ 
(Figure  2A)  the  impact  damage,  at  least  as  observed  on  the  surface, 
consists  of  a small,  well-formed  indentation.  Indentations  were  observed 
at  velocities  as  low  as  9 ms  At  71.4  ms  * (Figure  2B)  the  impact 
damage,  as  observed  on  the  surface,  consists  of  an  indentation  containing 
circumferential  cracks  and  surrounded  bv  radial  cracks.  At  still  higher 
velocities,  lateral  vent  cracks  extend  to  the  surface  leading  to  chipping 
between  the  radial  cracks. 

Observation  of  indentations  at  velocities  as  low  as  9 ms  \ means 
that  the  yield  stress  has  already  been  exceeded  at  this  velocity.  The 
stresses  were  estimated.  Static  loads  were  applied  to  the  surfaces 
using  the  same  spheres  as  indenters.  The  contact  radii  and  indentation 
radii  were  measured.  Then,  the  impact  loads  were  estimated  based  on  the 
assumption  that  the  impact  and  static  loads  of  the  same  magnitude  cause 
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Figure  2 Impact  Damage  in  As  Fired  TT-Zr02  Impacted  by 
1.588  mm.  Diameter  Tungsten  Carbide  Spheres  (82X) 


Indentations  of  the  same  size.  At  the  lowest  loads  at  which  indentations 
were  observed,  the  average  pressure  (p^)  on  the  indentation  was  about 
1 ClPa . Gilman^^  has  pointed  out  that  the  ratio  of  the  hardness  to  the 
yield  stress  increases  in  nonmetallic  crystals  depending  upon  the  relative 
difficulty  of  flow  on  the  various  glide  systems.  Therefore,  it  is  not 
possible  to  calculate  a definite  value  of  yield  stress  from  the  indenta- 
tion pressure.  According  to  the  Hertzian  stress  distribution,  the  radial 
stress  at  the  circumference  of  the  contact  area  is  0.166  p^  which  corresponds 
to  about  1.2  Gi’a  in  this  case.  This  localized  stress  is  about  three  times 
the  measured  flexural  strengths  of  this  material.  Therefore,  it  is  reason- 
able to  expect  flaw  growth  and  circumferential  crack  formation  under  the 
Influence  of  these  stresses. 

Energy  losses  during  impact 

For  impacts  of  tungsten  carbide  spheres  on  TT-ZrO.,,  the  coefficient 

of  restitution  decreases  approximately  linearly  with  impact  velocity  as 

shown  in  Figure  3.  The  slope  of  this  curve  contrasts  strikingly  with  the 

slope  obtained  at  room  temperature  for  impacts  of  glass  spheres  on  glass 

(12) 

plates  which  show  only  a very  slowly  decreasing  coefficient  of  restitu- 
tion. However,  the  slope  is  similar  to  that  observed  for  glass  at 
temperatures  above  the  softening  point  where  substantial  indentation 
occurs  because  of  viscous  flow^'\  This  similarity  is  evidence  that 
plastic  deformation  at  the  impact  site  is  responsible  for  much  of  the 
observed  decrease  in  the  coefficient  of  restitution.  The  aging  treatments 
have  little  effect  on  the  coefficient  of  restitution. 

The  kinetic  energies  of  the  impacting  and  rebounding  spheres  were 
calculated  anil  used  to  calculate  the  kinetic  energy  losses.  These  losses 


found  analytically 
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are  plotted  as  the  solid  line  In  Figure  4.  Hunter 
that  the  stress  wave  energy  induced  by  elastic  Impact  varied  as  the  2.6 
power  of  the  impact  velocity.  Evidently,  the  stress  waves  account  for 
all  or  almost  all  of  the  energy  losses  occurring  during  room  temperature 
impacts  of  glass  spheres  on  glass  surfaces^1 in  subsequent  experiments, 
it  was  found  that,  even  when  identations  were  formed  by  impacts  of  glass 
spheres  on  glass  surfaces  at  temperatures  above  the  softening  point,  the 
kinetic  energy  losses  still  varied  approximately  as  the  2.6  power  of  the 
impact  velocity . Therefore,  kinetic  energy  loss  was  plotted  vs.  the 
impact  velocity  on  a log- log  plot  in  Figure  5.  Clearly,  the  kinetic 
energy  loss  is  again  observed  to  vary  approximately  as  2.6  power  of  the 
impact  velocity  even  though  an  indentation  mechanism  contributes  substan- 
tially to  the  losses. 

Because  the  velocity  of  the  sphere  was  reduced  to  zero  by  impact  on 

a stationary  target,  conservation  of  energy  requires  that  the  incoming 

kinetic  energy  (U^)  is  equal  to  the  energy  transferred  to  the  target. 

Assuming  that  the  energy  transferred  to  the  target  consists  of  indentation 

energy  (II.)  and  stress  wave  energy  (U  ) the  energy  balance  can  be  written 
1 sw 


u 


k 


+ U 

sw 


(1) 


However,  consists  of  two  parts,  the  energy  of  plastic  deformation  (ll  ) 

and  the  energy  of  elastic  recovery  (UJ  which  accelerates  the  rebounding 

sphere.  The  kinetic  energy  loss  on  Impact  (AU,  ) is  11,  - l!  . Therefore, 

k k e 

subtracting  U from  both  sides  of  Equation  (1)  yields 


Energy  Loss- 
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Figure  4 Energy  Loss  vs.  Impact  Velocity  for  TT  Zr02  Impacted 
by  1.588mm.  Diameter  Tungsten  Carbide  Spheres 


I ■ 


Impact  Velocity -ms“' 


Figure  5 Kinetic  Energy  Loss  vs.  Impact  Velocity  for  TT  Zr02  Impacted 
by  1.588mm.  Diameter  Tungsten  Carbide  Spheres 
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AUk  = Hi  - Ue  + Usw  (2) 

By  comparing  estimates  of  AUj.  with  the  sum  of  the  terms  on  the  right  hand 
side  of  liquation  (2)  one  can  determine  whether  or  not  the  important  loss 
mechanisms  have  been  accounted  for. 

was  estimated  by  static  indentation  tests.  Tungsten  carbide  spheres 
were  loaded  to  successively  higher  loads  and  the  radius  of  the  indentation 
was  measured  at  each  step.  Assuming  a spherical  indentation,  the  penetra- 
tion (or  depth)  was  calculated  for  each  load  and  a load  vs.  penetration 
curve  was  plotted.  Integration  of  the  load  vs.  penetration  curve  yielded 
an  indentation  energy  vs.  penetration  curve. 

Usw  was  assumed  to  be  the  only  energy  loss  mechanism  at  velocities 
near  zero.  Usw  near  zero  velocity  was  estimated  by  extrapolating  the  curve 
in  Figure  J to  the  vertical  axis  and  calculating  the  kinetic  energy  loss. 

At  low  velocities  Usw  is  approximately  14%  of  the  impact  energy.  At  higher 
impact  velocities  Usw  was  assumed  to  be  a constant  fraction  of  the  kinetic 
energy  determined  by 

U = (1  - e“)  \ mv“  (J) 

SW  o Z 

in  which  e is  the  coefficient  of  restitution  near  zero  velocitv,  m is  the 
o 

mass  of  the  sphere  and  v is  the  velocity. 

The  results  of  calculations  for  individual  impacts  are  plotted  as 
open  circles  in  Figure  4 and  representative  values  are  given  in  Table  I. 

Even  though  the  results  are  scattered,  the  above  calculations  show  that  it 
is  reasonable  to  assume  that  deformation  and  stress  wave  energies  comprise 
the  major  energy  losses.  At  high  impact  velocities  the  energy  of  plastic 
deformation  (1)^)  becomes  the  dominant  loss  mechanism,  contributing  more 
than  80%  of  the  total. 
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TABLE  l 

Partitioning  of  Impact  Energy 
(As  Fired  TT-ZrO.,) 


Impact  Velocity  ms 

13.9 

71.4 

85.1 

Joules 

Impact  Energy  (U^) 

.00307 

.0811 

.1150 

Rebound  Energy 

.00225 

.0345 

.0162 

Energy  Loss  (ARE) 

.00082 

. 0466 

. 0998 

Indentation  Energy  d'j) 

.00167 

.0525 

.0974 

Plastic  Deformation  Energy  (U  ) 

P 

.00018 

.0180 

.0812 

Stress  Wave  Energy  (U  ) 

sw 

. 00044 

.0117 

.0166 

j 


Impact  damage  and  st  rengtli  dejjradat  ion 
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Goodier,  as  cited  bv  Evans  , found  that  the  plastic  penetration 
(indentation  depth)  q of  a projectile  can  be  expressed  for  small  penetra- 
t ions  by 


where  r is  the  radius  of  the  projectile,  v is  the  projectile  velocity, 

P is  the  projectile  density,  H is  the  hardness  of  the  target  and  6 is 
a constant.  In  the  present  case  the  projectile  and  target  properties 
are  constant  and  can  be  combined  with  to  yield  the  prediction  that  the 
penetration  is  proportional  to  the  velocity. 

The  penetration  depths  were  calculated  from  the  indentation  radii 
using  the  assumption  of  spherical  indentations  and  plotted  versus  the 
impact  velocity  as  shown  in  Figure  b.  The  penetration  increases  linearly 
with  velocity  as  expected  and  may  be  extrapolated  through  the  origin. 

This  observation  implies  that  the  1.588  mm  tungsten  carbide  spheres  act 
as  sharp  indenters,  i.e.  the  yield  stress  is  exceeded  at  all  impact  velo- 
cities . 

At  impact  velocities  over  about  25  ms  ^ radial  cracks  were  observed 
in  all  cases.  The  number  of  radial  cracks  was  counted  and  is  plotted 
versus  the  indentation  radius  in  Figure  7.  Extrapolation  of  the  sparse 
data  indicates  that  at  Indentation  radii  less  than  130-150  vim  one  should 
not  observe  radial  crack  formation.  In  fact,  indentations  with  radii  of 


107-124  pm  did  not  show  radial  crack  formation. 


As  Fired 
Under  Aged 


UUUJ-UOi;OJ49U9d 


o> 

iZ 


o 100  200  300  400 

Indent  Radius-pm 

Figure  7 No.  of  Radiol  Cracks  vs.  Indent  Radius  for  Impacts 
of  1.588  mm.  Diameter  Spheres  on  TT  Zr02 
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hunge  and  Kvatis^  Impact  Oil  yttrl.i  stabilized  zlrconia  using  400 

pm  tungsten  carbide  spheres  at  velocities  up  to  750  ms  1 and  observed  no 
radial  crack  tormatlon.  The  absence  ot  radial  cracks  mav  possibly  be 
explained  bv  Insul I Icient  strain  at  the  bi>undary  lor  verv  small  Indenta- 
t Ions . 

The  number  ol  radial  cracks  seems  to  decrease  with  Increasing  aging 
time.  This  decrease  may  have  an  Important  el  led  on  the  strength  degrada- 
tion as  will  be  discussed  later. 

The  Length  ol  the  longest  radial  crack  In  each  specimen  was  measured. 
The  radial  crack  lengths  are  plotted  versus  the  impact  velocities  In 
Figure  8.  The  radial  crack  lengths  ot  specimens  with  the  various  aging 
times  are  nearly  equal  at  given  velocities  up  to  about  70  ms  *.  This 
observation  Is  oont rary  to  expectations  based  on  the  Kj  values  ot  the 
materials  and  will  be  discussed  In  more  detail  later. 

The  remaining  strengths  of  l lie  Impacted  specimens  were  measured  (one 
bar  for  each  test  condition)  by  three  point  Loading  on  a one  Inch  span  with 
the  Impact  site  on  the  tensile  surface  under  the  center  Load  point.  The 
strengths  decreased  with  increasing  impact  velocity,  as  expect  oil  (Figure  '>) . 
The  shapes  ot  the  curves  for  the  various  aging  treatments  were  similar. 

Fach  treatment  maintained  its  relative  strength  advantage,  approximately, 
at  each  velocity. 

Kj  , values  were  estimated  from  the  fracture  stresses  and  measured 
t law  sizes,  for  a small  number  of  specimens  with  rather  small  flaw  sizes 
wlLlt  the  results  shown  In  Table  I.  The  Kj  , values  Increase  with  aging 
time  as  expected  but  the  values  are  somewhat  lower  than  those  obtained 


|j 

H 

i 

i 


by  the  SKN11  method.  In  Interpreting  these  results  it  should  be  remembered 
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Figure  9 Flexural  Strength  vs.  Impact  Velocity  for  TT  Zr02  Impacted 
by  1.588mm.  diam.  Tungsten  Carbide  Spheres 
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that  these  values  are  based  on  preexisting  flaw  measurements  but 
substantial  subcritical  crack  growth  may  have  occurred  leading  to  critical 
flaws  larger  than  those  measured.  Therefore,  these  estimates  of  K 

X c 

may  be  smaller  than  the  actual  values.  To  obtain  final  estimates  of  K , 

L L 

the  SENB  results  which  are  thought  to  be  too  high  and  those  based  on  the 
flaw  measurements  which  are  thought  to  be  too  low  were  averaged  yielding 
the  results  which  are  included  in  Table  I. 

As  noted  in  Table  I,  Gr if f ith- Irwin  fracture  mechanics  can  be  used 
to  calculate  the  fracture  stresses  of  specimens  with  flaws  of  known  sizes 
and  shapes.  The  basic  equation  can  be  written  as 


KIC=Y  °FC 


1/2 


(5) 


in  which  Y is  a factor  which  accounts  for  the  flaw  shape  and  proximity 

to  the  surface,  0 is  the  fracture  stress,  and  C is  the  flaw  size, 

F 


usually  the  flaw  depth.  Lawn  and  Marshall 


(19) 


have  related  the  depth 


of  damage  induced  by  contact  stresses  to  the  load  (P)  and  K by 

lv 

dimensional  analyses  as 


K. 


IC 


(6) 


in  which  X is  a dimensionless  contact  constant.  Combining  these  equations 
to  eliminate  C yields 


K. 


4/3 


a _ IC 

°f  - TTT773 


(7) 
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Equation  (7)  leads  one  to  expect  that , if  P is  constant  for  the  various 
materials  impacted  at  a particular  velocity,  0 « K Using  the  average 

values  in  Table  I,  one  would  expect  the  fracture  stresses  of  the  peak- 
aged  specimens  to  be  slightly  more  than  twice  those  of  the  as-fired 
specimens  but  the  differences  are  less  than  that.  This  difference  can 
be  accounted  for  by  two  factors.  First,  the  loads  exerted  on  the  under- 
aged and  peak-aged  materials  at  a given  velocity  are  slightly  greater 
(15-25%)  than  the  loads  exerted  on  the  as-fired  material.  Second,  the 
smaller  number  of  radial  cracks  in  the  peak-aged  material  as  compared  to 

the  as-fired  leads  to  higher  stress  intensity  values  at  a given  radial  crack 

* 

length  so  that  the  radial  cracks  in  the  peak-aged  material  will  tend  to 
be  longer  than  otherwise  expected.  The  effect  of  the  number  of  radial 
cracks  on  the  stress  intensity  factor  is  not  accounted  for  in  the  theory. 

Based  on  Equation  (7)  one  would  expect  the  remaining  strength  of  a 
particular  material  impacted  by  a particular  projectile  to  vary  linearly 
with  P In  the  present  case,  one  would  expect  three  curves  with  dif- 

ferent slopes  because  the  three  materials  have  different  K values. 

The  sparse  data  (Figure  10)  are  evidently  not  sufficient  to  demonstrate 
this  last  effect,  but  the  data  taken  together  show  a linear  relationship 
although  the  curve  does  not  pass  through  the  origin.  Using  reasonable 

estimates  of  Y and  x»  O’  * 1.11,  \ = .102),  and  the  measured  slope  yields 

1/2 

an  estimate  of  K^,  which  was  8.6  MPam  , a value  in  the  range  of  the 
measured  values. 


•k 

It  is  well  known  that  the  presence  of  other  cracks  can,  in  some  cases, 
reduce  the  stress  intensity  factor  at  the  tip  of  a particular  crack. 

The  authors  are  not  aware  of  a solution  for  radiating  semi-circular  cracks. 
However,  Tada,  Paris  and  lrwin(20)  contains  a solution  for  radiating  straight- 
through  cracks.  In  this  case,  for  more  than  ten  radial  cracks  the  ratio  of 
the  stress  intensity  factor  for  n radial  cracks  compared  with  that  for  one 
such  crack  is  approximately  2/  / n. 


Strength  After  Impoct  for  TT  ZrO 
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IV.  Conclusions 

IT  ZrO , is  resistant  to  localized  impact  damage  bv  glass  and  steel 
spheres.  I mm  diameter  glass  spheres  impacting  at  velocities  up  to 
124  ms  * and  1.59  nun  diameter  steel  spheres  impacting  at  velocities  up 
to  198  ms  did  not  cause  visible  damage.  Harder  tungsten  carbide 
spheres,  1.59  nun  diameter,  indented  the  surfaces  at  impact  velocities 
as  low  as  9 ms  showing  an  elastic  plastic  response.  At  higher  impact 
velocities  radial,  circumferential  and  lateral  vent  cracks  were  observed, 
in  addition  to  the  indentations. 

The  energy  loss  during  impact  was  analyzed.  The  work  of  plastic 
deformation  accounts  for  most  ot  the  energy  loss  at  high  velocities.  Most 
of  the  remaining  energy  loss  can  be  accounted  for  by  the  stress  wave  energy. 

The  remaining  strengths  of  the  TT  ZrO,  specimens  with  three  aging 
treatments  decreased  with  Increasing  impact  velocity,  as  expected.  At 
a particular  impact  velocity,  the  remaining  strength  increased  with 
increasing  aging  time  fK^,)  hut  the  increase  was  somewhat  less  than  expected. 
This  difference  can  be  accounted  for  by  the  slightly  greater  loads  on 
the  higher  K specimens  and  the  formation  ot  fewer  radial  cracks.  Higher 
loads  and  reduced  number  of  radial  cracks  yield  higher  stress  intensity 
factors  and  greater  crack  propagation  than  otherwise  expected. 
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Abstract 


Crack  propagation  and  branching  were  investigated  in  calcia  partially 
stabilized  zirconia  aged  at  1300°C  for  various  times.  The  stress  intensity 
factor  at  crack  branching  increased  with  increasing  aging  time.  The  strain 
intensity  criterion  for  crack  branching  was  used  to  estimate  the  effective 
crack  tip  radius  at  branching  and  the  increase  in  this  radius  attributable 
to  the  transformation  toughening  mechanism. 
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1.  Introduction 

In  brittle  materials  there  is  an  orderly  progression  of  fracture  features 
along  lines  drawn  radially  in  the  fracture  surface  from  the  fracture  origin. 

These  features  include  a flat  region  usually  referred  to  as  the  mirror,  a 
mist  region,  a region  of  radiating  ridges  and  valleys  called  hackle  and  a 
boundary  at  which  large  scale  crack  branching  occurs.  In  this  paper  large 
scale  crack  branching  is  simply  called  crack  branching  and  is  to  be  distinguished 
from  microbranching.  Crack  branching  occurs  at  a particular  value  of  stress 
intensity  factor  and  is  a response  to  large  scale  stresses  and  macro  properties. 
Micro-branching  is  considered  to  be  a response  to  localized  stresses  and 
material  property  variations. 

The  fracture  stress  (o  ) and  the  radius  to  the  crack  branching  boundary 

r 

(1-4) 

(C„)  can  be  related  bv 
n 

kb  ■ yVb1/"’  <» 

in  which  is  the  stress  intensity  factor  at  crack  branching  and  V is  the 

appropriate  geometrical  factor  for  the  crack  at  branching.  IT  can  be  determined 

B 

-1/2 

bv  plotting  a vs  C “ and  determining  the  slope  which  is  K ,/Y. 

b B b 

The  fracture  toughness  o'"  zirconia  can  be  increased  by  proper  control 
of  composition  and  heat  treatment ^ When  calcia  partially  stabilized 

zirconia  is  aged  at  intermediate  temperatures  ( 1200-1500°C)  in  the  cubic- 
tetragonal  phase  field,  the  tetragonal  precipitate  in  the  much  larger  cubic- 
crystals  is  coarsened^ . At  lower  temperatures  the  tetragonal  phase  becomes 


increasingly  susceptible  to  transformation  to  the  monoclinic  phase.  Tensile 
stress  favors  this  transformation  because  the  monoclinic  phase  is  less  dense 
than  the  tetragonal.  Therefore,  when  cracks  propagate  through  the  material 


tlu1  tensile  stress,  concentrated  at  the  crack  tip,  causes  the  tetragonal 

(9 ) 

crystals  within  a tow  micrometers  ot  tlu*  crack  to  transform  to  monoclinie 
This  transformation  increases  t ho  crit  leal  st  less  Inlonsltv  factor  (Kj(1)  so 
that  tlu>  malorlal  is  some!  imos  rail  oil  l ranst  ormat  1 on  t oughenod  zircon  in 
( TTZrl)  ,) . 


Tho  mechanism  rosponsihlo  lor  tho  Inoroaso  in  fracture  toughness  will) 
increasing  transformation  from  tho  tetragonal  t v>  tho  monoclinic  is  not  com- 
pletely  understood  * ^ \ In  some  cases  tho  Individual  crystals  are  con- 
sidered ici  act  as  harriers  t o crack  propagation^*  ^ or  as  a source  ot  micro- 
crack  formation*'  ’ ' and,  in  other  cases.  It  is  assumed  that  the  expansion 


that  occurs  with  format  ion  ol  t lie  monocllnlc  phase  absorbs  energy  and  acts 


to  reduce  elastic  strain  at  the  crack  tip 


ll«>,17) 


The  lactors  determining  K in  a particular  ceramic  material  and  the 
relationship  ol  to  K also  are  not  understood  completely.  In  most  ceramics, 
the  K does  not  vary  much  Iroin  the  value  one  might  expect  based  on  the 
variation  ol  Young's  modulus.  This  lack  ol  variability  has  made  It  difficult 
to  study  toughening  mechanisms  and  the  relationship  ol  to  K in  ceramics, 
because  ol  the  high  strength  and  I incline  toughness  ot  TT7.rO.,  and  the  tact 
that  these  properties  can  be  varied  bv  heal  treatment,  with  little  change  in 
Young's  modulus,  it  was  ol  interest  to  studv  crack  propagation  and  crack 
branching  in  TT-ZrO., . 


St) 


I I . Proeedures 

Calola  part  ial  ly  stahi 1 Ued  ZrO.,  with  four  l .100*1;  aging  t reatmeuts 

(as-tired,  0 hours;  under-aged , 10  hours;  peak-aged , 48  hours;  and  over-aged, 

*>/  hours)  was  propon'd  at  I'SIKO  in  t ho  form  ol  root  angular  bars,  1 x 10  \ 40mm, 

with  ground  (‘>0  grit)  surlaoes.  At  I'eramio  Finishing  Company  ono  largo  surtaoo 

ot  uuih  bai  was  polishod  on  a papoi  oovoroil  rotating  lap.  Polishing  romovod 

most  ot  t ho  monoellnlo  phaso  induood  in  t ho  surtaoo  by  grinding  and  tho 

(19) 

rosulting  oomprossivo  surtaoo  lavoi  . Diamond  polishing  was  moro  effective 

than  alumina  polishing  in  reducing  tho  monoolinio  oontont  ot  tho  surtaoo. 

l'he  specimens  analyzed  in  this  invost igat  ion  woro  originally  traoturod 

(19) 

In  an  invost  igat  ion  ot  strongth  ilogr.ulat  ion  duo  to  localized  impaot  damage 
I'ho  spooimons  that  woro  analyzed  woro  those  that  traoturod  at  relatively  high 
stresses;  that  is,  tho  oontrols  that  woro  not  impaotod  and  tho  spooimons  with 
minimal  impaot  d.onago.  Only  those  spooimons  showed  oraok  hranohing.  the 
spooimons  woro  tested  In  t lexure  by  throe  point  loading  on  a ono  iueh  span. 

I’ho  traoturo  surtaoos  ot  tho  spooimons  wore  oharao  t or  i zed  . I’ho  tollowing 
eritorla  woro  used  to  looato  tho  traoturo  origins: 

l.  tin'  eentei  ot  area  ot  rot  looting  spots  wtiioh  woro  observed  bv 
optioal  miorosoopy  with  tho  spooimon  illuminated  at  tho  propot 
nng I o . 

fin'  venter  ot  somi-eirelo  termed  by  flakes  of  unitorm  size. 

1.  The  int orsoo t ton  ot  linos  extended  through  the  mirror  region 
para  1 t e l to  t he  hao k 1 o . 

4.  fho  ooutor  ot  tin'  smooth,  relatively  flat  area  (traoturo  mirror), 
file  tollowing  features  were  oharao t or  l zed ; traoturo  origin,  ret  looting  spots, 
oraok  hranohing  and  mlorobranoh ing.  The  radii  at  oraok  hranohing  wore  measured 
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using  an  optical  microscope  with  a calibrated  eyepiece.  The  measurements  were 
made  on  the  tensile  surface  where  the  crack  propagates  in  a uniform  stress 
t le Id. 

K was  calculated  as  described  in  Section  1.  Y was  assumed  to  be  1 . !18 

which  is  the  value  determined  by  Smith,  Emery  and  Kobayashi  for  the 

greatest  Kj  at  the  boundary  of  a semi-circular  surface  flaw.  In  some  cases, 

the  fractures  originated  near  corners.  One  might  question,  whether  or  not 

the  presence  of  the  additional  free  surface  would  significantly  increase  the 

( 21) 

stress  intensity  factors  in  these  cases.  However,  as  shown  by  Liu'  , this 
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111.  Results 

Characterization  of  fracture  surfaces 

The  reflecting  spots  were  the  principal  features  used  to  locate  the 

fracture  origins.  In  many  cases,  the  locations  of  the  origins  were  confirmed 

using  the  other  criteria  listed  in  the  previous  section.  The  reflecting  spots 

are  cleavage  faces  ol  the  zirconia  grains  and  are  similar  to  those  observed  at 

(22,23) 

fracture  origins  in  alumina  ceramics  . However,  the  grain  size  of  the 

zirconia  is  very  large  (-  50pm)  so  the  reflecting  spots  are  correspondingly 
large.  In  some  cases  the  ret  laeting  surfaces  appear  to  consist  ot  fan  shaped 
segments  that  are  oriented  at  slightly  varying  angles  to  the  incident  light 
so  that  slight  changes  in  spec ime  i orientation  are  necessary  to  bring  the 
individual  segments  into  reflecting  position.  Because  the  grain  size  is  so 
large,  only  a few  reflecting  spots  are  observed  at  the  fracture  origin. 

In  a small  number  of  cases,  the  flaw  at  the  fracture  origin  was  large 
enough  to  be  identified  and  characterized.  In  most  cases,  the  particular 
flaw  at  which  the  fracture  originated  was  not  identified. 

The  fracture  surface  near  a fracture  origin  is  illustrated  in  Figure  1. 
The  fracture  origin  was  a large  pore.  This  fracture  surface  was  chosen  because 
identification  of  the  flaw  at  the  fracture  origin  was  certain.  The  fracture 
is  almost  entirely  t ransgranular . The  fracture  origin  is  surrounded  by  a 
relatively  flat  region  containing  some  radiating  lines  that  may  be  steps 
formed  by  twisting  of  tin'  crack.  The  flat  region  is  surrounded  by  a wnvv 
region  beginning  at  about  0.4mm  from  the  fracture  origin.  The  wavy  region 
becomes  increasingly  rough  as  the  distance  from  the  fracture  origin  increases. 

The  interiors  of  pores  and  grain  boundaries  exposed  by  fracture 
frequently  show  verv  uniform  surface  roughness  in  limited  areas  (Figure  2). 
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This  roughness  is  believed  to  result  from  the  presence  of  the  precipitate  in 
the  larger  cubic  grains.  The  presence  ol  the  grain  boundary  roughness  may 
explain  the  resistance  to  intergranular  fracture  tint  results  in  the  almost 
complete  absence  of  intergranular  fracture  near  fracture  origins. 

At  distances  ranging  from  about  one  to  tour  millimeters  from  the  fracture 
origin,  depending  on  fracture  stress  and  aging  time,  the  first  crack  branching 
was  observed.  Unlike  most  ceramics  in  which  branching  is  observed,  the 
TT-ZrO,  specimens  fractured  in  two  halves  rather  than  forming  a number  of 
wedge  shaped  fragments  in  addition  to  the  two  end  pieces.  Apparently  the 
increased  microbranching  and  at  crack  arrest  in  the  TT-ZrO.,  prevent  the 
branches  from  propagating  to  the  point  that  the  wedge  shaped  pieces  can  separate 
from  the  two  halves.  The  branching  can  be  observed  in  the  tensile  surfaces  of 
the  specimens  under  inclined  illumination.  The  edges  of  the  cracks  appear  to 
be  slightly  raised  so  that  they  can  be  seen  clearly.  Frequently,  two  or  more 
branches  are  observed  at  nearly  equal  intervals  from  the  fracture  origin. 

The  branches  were  also  observable  on  the  fracture  surface.  Because 
the  TT-ZrO,,  is  rather  transparent,  light  is  transmitted  through  the  edges  of 
the  wedges  formed  by  branching  and  reflected  from  the  far  surface  so  that  the 
branches,  near  their  origins,  appear  as  irregular  light  colored  bands. 

A branch  formed  in  the  tensile  surface  of  a peak-aged  specimen  fractured 
at  a relatively  high  stress  of  51bMPa  is  illustrated  in  Figure  3.  In  the 
photograph  the  edge  formed  bv  the  main  crack  is  approximately  horizontal  and 
the  branch  slopes  downward  at  a branching  angle  of  about  25°.  The  wedge  shaped 
piece  formed  by  the  branch  did  not  separate  from  the  rest  of  the  specimen.  The 
microbranching  along  the  branch  is  evident. 
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In  the  optical  microscope,  raised  edges  of  the  cracks  in  TT-ZrO.^  cast 
a shadow  that  makes  the  cracks  visible.  Various  attempts  have  been  made  to 
observe  these  features  in  the  SEM  but  little  progress  has  been  made.  A typical 
result  is  illustrated  in  Figure  4 which  shows  cracks  that  propagated  across  a 
relatively  smooth  surface  with  microbranching  at  the  centers  of  the  photographs. 
At  the  higher  magnification  and  45°  incidence  it  appears  that  the  material  in 
the  wedge  between  the  microbranches  is  raised  slightly  above  the  original 
surface.  A more  ridge-like  feature  extending  along  the  cracks  caused  by  the 
volume  expansion  during  the  phase  transformation  had  been  anticipated  but  this 
has  not  yet  been  observed. 

Microbranching  is  common  in  these  materials.  As  the  crack  grows 
microbranching  is  first  observed  near  the  first  large  scale  branching.  There- 
after, microbranching  is  observed  more  or  less  continually  along  the  branches. 
These  observations  suggest  that  the  stress  intensity  factor  for  microbranching 

is  close  to  K . Furthermore,  the  absence  of  microbranching  near  fracture 
n 

origins  indicates  that  other  mechanisms  are  responsible  for  the  toughening. 

It  should  be  emphasized  that,  although  fracture  mechanisms  that  occur  when 

> K1C  are  interest ^n8  i'1  their  own  right,  they  cannot  provide  evidence  about 
toughening  mechanisms.  This  evidence  can  only  be  obtained  by  observations  of 
the  fracture  surface  where  KT  = K.„. 

Crack  branching  measurements  and  calculation  of  Kk 

The  crack  branching  radii  were  measured  and  plotted  as  suggested  by 

equa.  (1)  in  Figure  5.  The  data  for  the  various  aging  times  fall  into  three 

groups  yielding  slopes  that  increase  with  aging  time.  K was  calculated  from 

B 

the  data.  The  average  K values  for  various  aging  times  and  surface  finishes 

n 

are  given  in  Table  I. 


Figure  5 Fracture  Stress  vs.  cB“l/2  for  TT-Zr02  Specimens  with 
Various  Aging  Treatments  (Combined  Data  for  Three 
Different  Surface  Finishes). 


90 


cx 

3 

O 

u 

ae 

u 

0) 

a. 

c 

< V 
B 

•H 

u 

0) 

ex 

CO 

OJ 

c 

o 


■K 

•K 


91 


IV.  Discussion  of  Results 


Evaluation  of  the  effect  of  residual  compressive  surface 
stresses  on  crack  branching 

Analysis  of  crack  branching  data  in  cases  in  which  large  scale  residual 
stresses  can  be  present  must  begin  with  a determination  of  the  importance  of 

(3 

these  stresses  and  a decision  on  how  to  account  for  this  added  complication 

The  presence  of  residual  compressive  surface  stresses  (0  ) decreases  the 

nominal  applied  stress  (o  ) acting  at  the  fracture.  Therefore,  the  fracture 

stress  (o  ) is  algebraic  sum  of  the  stresses  as  in 
F 


+ cr_ 


(2) 


where  compressive  stresses  are  negative.  Substituting  in  equa.  (1)  yields 


kb  ■ V<<,N 


0 )C 

r'  b 


1/2 


(3) 


Flexural  strength  measurements  consistently  show  that  the  strengths  of 

TT-ZrO^  specimens  increase  when  they  are  ground  and  decrease  when  they  are 

polished This  same  effect  was  observed  in  the  present  specimens  as  shown 

in  Table  II.  The  greatest  strength  difference  (98MPa)  is  observed  for  the 

peak  aged  specimens  which  are  most  susceptible  to  the  tetragonal  to  monoclinic 

(23) 

phase  transformation.  Pasco  and  Garvie  attribute  the  strengthening  to  a 
compressive  surface  layer,  20-30ym  in  depth,  that  is  generated  by  abrasion  and 
removed  by  polishing. 

The  K data  in  Table  I which  were  calculated  using  c„  show  little  variation 

D N 

with  the  surface  finish.  This  result  was  surprising  in  view  of  the  rather  large 
strength  differences  but  it  may  be  caused  by  the  fact  that  the  stressed  layers 
are  very  thin.  In  any  case,  for  the  remainder  of  the  analysis  it  is  assumed  that 


values  for  the  TT-ZrO^  materials  were  estimated  previously  by 

averaging  results  of  single  edge  notched  beam  measurements  and  results  calculated 

(19) 

from  flaw  sizes  and  fracture  stresses  . These  K values  and  the  present  K 

XC  B 

values  are  compared  in  Table  III.  Probably,  the  K values  were  not  determined 

i-  Li 

accurately  enough  to  merit  detailed  comparison,  as  also  seems  to  be  the  case 

/ r\  I nr  \ 

even  for  intensively  investigated  ceramics  such  as  H.  P.  silicon  nitride'  ’ 


Nevertheless,  it  is  clear  that  increases  as  K increases.  The  Young's 
moduli  of  the  materials  vary  only  slightly  with  aging.  Because  K = (2EY  ) ‘ 

J.  L*  C 

where  Yc  is  the  critical  fracture  energy,  it  is  clear  that  the  increase  of 

with  aging  is  the  result  of  increased  fracture  energy. 

(4) 

The  strain  intensity  criterion  for  crack  branching  can  be  used  to 

estimate  the  effective  radius  (r  ) °f  the  tip  of  the  propagating  crack  at 

branching.  To  do  this  the  strain  intensity  at  crack  branching  (K  /E)  is 

B 

calculated  and  then  squared  to  obtain  the  estimate  of  r . Earlier  results 

e 

for  a large  number  of  brittle  ceramics  yielded  an  average  r^  of  11  x 10  ^m,  a 

value  about  10  times  atomic  dimensions.  The  present  estimates  for  TT-ZrO^ 

range  from  8 to  20  times  as  great  as  those  previously  calculated  for  more 

conventional  brittle  ceramics. 

(18) 

Lange  investigated  crack  propagation  in  100%  cubic  zirconia  in  which 

no  tetragonal  to  monoclinic  transformation  occurs  and  observed  a value  of 

1/2 

3MPam  . Comparison  of  this  value  with  the  value  for  the  as-fired  material 

. 1/2 

in  this  investigation  (K^  = 5.1MPam  ) indicates  that  even  in  the  as-fired 
material  transformation  toughening  is  important.  If  one  assumes  that,  in  the 


absence  of  the  transformation  mechanism  r^  of  ZrO^  is  the  same  as  the  average 
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B — 5 1/2 

value  tor  the  brittle  ceramics  investigated  previously  (—  = 3.3*10  m so 

h 

that  r = 10.9*10  ^m) , one  can  estimate  the  increase  in  r due  to  the 
e e 

-9  -9 

transformation  mechanism  (tor  example,  for  peak-aged  TT-ZrO^  21.7*10  -1.05*10 

-9 

20.7*10  m) . It  is  an  interesting  question  whether  or  not  this  change  in  r^ 
could  be  estimated  from  the  volume  expansion  of  the  tetragonal  to  monoclinic 
transformation,  the  fraction  of  tetragonal  phase  originally  present  and  the 
variation  of  the  percentage  of  transformation  occurring  with  distance  from  the 
crack  surface.  This  estimate  lias  not  been  attempted  here  because  of  insufficient 
information  on  the  last  two  items. 

One  might  question  why  it  is  desirable  to  invoke  the  strain  intensity 

criterion  for  crack  branching  when  it  is  clear  from  available  evidence  such  as 

Figure  l of  reference  (13)  that  cracks  in  TT-ZrO.,  are  deflected  in  many  cases 

by  the  precipitate  particles,  indicating  that  the  particles  resist  crack 

(22) 

propagation.  The  point  is  that,  as  shown  recently  for  96%  alumina  , the 
increase  in  K^,  as  the  crack  grows  beyond  the  critical  flaw  boundary,  must  be 
accompanied  by  absorption  of  energy  by  a variety  of  nut  yet  identified  mechanisms 
The  crack  can  only  absorb  a given  amount  of  energy  in  going  around  a particular 
precipitate  particle.  The  remainder  must  be  absorbed  by  other  mechanisms 
perhaps  involving  extension  of  the  thickness  of  the  transformed  layer.  Thus, 
r_  continues  to  increase  as  the  crack  propagates. 

It  is  important  to  realize  the  r does  not  represent  the  actual  radius 
of  the  crack  tip  which  remains  small.  Neither  does  it  represent  the  radius 
of  the  zone  in  which  transformation  occurs.  Instead  it  represents  the  radius 
the  crack  tip  would  have  to  have  if  the  expansion  due  to  the  phase  transformation 
had  not  occurred  but  the  fracture  stress  and  fracture  energy  remained  the  same. 
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It  one  views  crack  branching  as  a fracture  process  that  occurs  when  the 
total  elastic  strain  at  the  crack  tip  reaches  a critical  value  sufficient  to 
support  propagation  of  two  cracks.  It  is  clear  that  mechanisms  such  as  phase 
transformations  occurring  with  an  Increase  in  volume,  and  plastic  deformat ion, 
will  act  to  reduce  the  total  elastic  strain  at  each  particular  crack  length, 
postponing  crack  branching  until  the  elastic  strain  intensity  increases  to 
the  point  that  crack  branching  can  occur. 
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